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FOREWORD 


This  report  was  prepared  by  Analytical  Services  &  Materials,  Inc.,  Hampton  Virginia  for 
WL/FIBEC,  Wright-Patterson  Air  Force  Base,  Ohio  under  contract  F33615-94-D-3212, 
“Structural  Integrity  Analysis  cind  Verification  for  Aircraft  Structures.”  The  contract 
program  manager  was  1st  Lt  Dave  Conley,  WL/FIBE.  The  government  project  engineer 
was  James  A.  Harter.  The  period  of  performance  for  this  report  was  1  Jan  95  through  30 
Sept  95. 


The  work  was  performed  under  project  (Delivery  Order  0003),  by  Analytical  Services  & 
Materials,  Inc.  personnel  located  at  the  WL/FIBEC  Fatigue  &  Fracture  Test  Facility, 
Bldg.  65,  Area  B,  Wright-Patterson  AFB,  OH.  This  effort  was  a  follow  up  to  previous 
work  performed  under  Delivery  Order  0001  [1].  The  Principal  Investigator  of  this 
research  was  Mr.  Kevin  L.  Boyd.  The  authors  of  this  report  were  Mr.  John  H.  Eisner,  Mr. 
Kevin  L.  Boyd,  and  Mr.  Srinivas  Krishnan.  Technical  inputs  were  submitted  by  Mr. 
James  A.  Harter  and  Mr.  Daniel  A.  Jansen. 


1. 


EXECUTIVE  SUMMARY 


This  research  included  baseline  crack  growth  data  for  Ti-6-4  and  4340-195ksi  steel, 
verification  testing  of  the  Ti-6-4  data,  testing  various  geometries,  and  AFGROW 
development. 

Baseline  test  data  were  produced  for  Ti-6-4  and  4340-1 95ksi.  Center  cracked  specimens 
were  tested  at  four  stress  ratios.  Load  shedding  enabled  crack  growth  rates  as  low  as 
10.0'®  to  be  recorded  for  both  materials.  Determining  crack  growth  rates  at  this  level  is 
important  for  life  predictions  since  much  of  the  life  is  accounted  for  when  cracks  are 
small. 

Open  hole  Ti-6-4  specimens  were  tested  to  verify  the  crack  growth  rate  data  generated 
and  to  test  AFGROW’s  predictive  abilities.  For  the  six  open  hole  specimens,  AFGROW 
had  2in  average  error  of  12%. 

Joint  specimens  were  also  tested  in  order  to  verify  AFGROW’s  abilities.  These 
specimens  were  modeled  with  finite  element  analysis  to  determine  the  load  transfer 
around  the  rivet  holes.  Once  this  was  accomplished,  predictions  were  run  and  compared 
to  test  data.  The  results  of  the  single,  multiple  in-line  joint,  and  multiple  staggered  joint 
specimens  were  good. 
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2. 


INTRODUCTION 


The  purpose  of  this  research  is  to  determine  fatigue  crack  growth  rates  for  Ti-6-4  and 
4340-195ksi  steel  and  also  to  verify  existing  crack  growth  databases  in  AFGROW  [2]  by 
testing  various  geometries  and  comparing  the  predictions  to  test  results.  AFGROW  is  a 
fatigue  crack  growth  prediction  code  developed  by  AS&M  and  the  Air  Force.  In  order  to 
make  AFGROW  more  useful,  common  structural  materials  need  to  be  included  in  the 
material  database  that  is  included  with  the  program.  The  materials  database  in  AFGROW 
is  based  on  simple  middle  tension  (MT)  or  compact  tension  (CT)  specimens.  The  majority 
of  data  are  generated  with  MT  specimens  where  data  are  sampled  for  as  wide  of  a  range 
of  crack  growth  rate  (da/dN)  as  possible.  Normally,  data  in  the  range  of  2.0  E-09  to  1.0 
E-03  in/cycle  are  sampled  using  a  decreasing  -  increasing  delta  K  method.  Of  course, 
great  care  is  taken  to  avoid  overload  effects  as  far  as  possible  and  data  taken  when  delta  K 
is  increasing  are  compared  to  that  taken  during  the  initial  decreasing  period.  The 
predictions,  which  rely  on  AFGROW’ s  database,  are  then  verified  with  test  data  that  were 
generated  with  different  and  more  complicated  geometries. 
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3. 


FINITE  ELEMENT  ANALYSIS 


The  three  bolted-joint  configurations  examined  in  this  report  were  analyzed  to  estimate  the 
amount  of  load  transfer  present  in  the  specimens.  This  information  was  required  so  that 
the  life  predictions  could  be  made  for  the  bolted-joint  configurations  using  the  AFGROW 
crack  growth  life  prediction  code  [2], 

For  this  study,  typical  double  shear  joint  specimens  were  designed  to  address  an  additional 
order  of  complexity  in  the  life  prediction  process,  or  cracks  growing  from  loaded  holes. 
These  “typical”  joint  configurations  were  designed  with  acceptable,  industry-standard  joint 
design  parameters  such  as  fastener  size  (0.25  in.  typ.),  pitch  distance  (1.0  in.  typ.)  and 
edge  distances  (0.75  in.  typ.).  The  single  fastener  test  specimens  were  16.0  inches  in 
length  and  1.0  inch  in  width.  The  plate  thicknesses  for  the  single  fastener  specimens  were 
0.25  inch.  The  multiple  fastener,  in-line  and  staggered  specimens  were  16.0  inches  in 
length  and  3.95  inches  in  width.  The  plate  thicknesses  for  the  multiple  fastener  specimens 
were  0.09  inch.  These  specimen  configurations  and  their  crack  locations  can  be  seen 
below  in  Figure  1 . 


Crack  Location 


Crack  Location  - 

Crack  Location  - 

V 

L 

•  •  • 

•  •  • 

;  • 

•  •  • 

•  #  • 

•  •  • 

•  • 

Single  Fastener 


Multi  Fastener 
(In-Line) 


Multi  Fastener 
(Staggered) 


Figure  1:  Bolted- Joint  Configurations 
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All  of  the  joint  configurations  modeled  in  this  study  were  of  the  symmetric  lap  splice,  or 
double  shear  configuration.  These  configurations  were  chosen  to  reduce  the  amount  of 
bending  present  in  the  specimens  during  testing.  During  the  testing  phase  of  this  study,  all 
specimens  were  fatigued  under  uniaxial  loading  conditions.  Also,  all  fasteners  were 
protruding  head,  small  clearance,  or  “neat”  fit  A286  steel  fasteners  that  were  lightly  finger 
tightened,  to  remove  the  effect  of  load  transfer  through  friction  or  “clamp-up”  and/or  pre¬ 
stress  effects  of  the  bolts. 

The  single  fastener  joint  was  assumed  to  have  100%  load  transfer  for  crack  growth  life 
prediction  purposes.  This  joint  configuration  was  verified  with  that  FRANC2D/L  finite 
element  code  [3]  and  is  shown  in  Figure  3.  Symmetry  conditions  were  assumed  both 
through  the  thickness  of  the  joint  and  along  the  specimen’s  length.  The  planes  of 
symmetry  assumed  in  the  analysis  can  be  seen  below  in  Figure  2. 


"•1 


Plane  of  Syinmetry 


Plane  of  Symmetry 


Figure  2:  Symmetry  Planes  Assumed  for  2D  Finite  Element  Analysis 


Because  FRANC2D/L  is  a  two  dimensional  layered  code  only  one  layer’s  mesh  can  be 
seen  at  a  time.  The  mesh  for  the  “joined”  member  is  shown  in  Figure  3.  The  “outlined” 
splice  plates  mesh  must  exactly  mirror  the  mesh  that  is  displayed  within  the  confines  of  the 
displayed  box.  The  mesh  of  the  splice  plate  is  not  shown  due  to  the  redundancy  of  the 
mesh.  The  material  properties  used  in  the  analysis  are  also  shown  in  Figure  3.  A  unit 
remote  load  was  applied  to  the  specimen  and  the  rivet  load  was  calculated  to  be  one  (1.0), 
signifying  100%  load  transfer. 
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Location  of  Rivet  Element 


Material:  7075-T7351  Aluminum 
E  =  10.3  E6  psi/in^ 
v  =  0.33 

Member  Thickness  =  0.125  in 
Splice  Plate  Thickness  =  0.25  in. 
Rivet  Stiflhess  =  6.06  xlO^  Ib/in 


Figure  3:  Single  Fastener  Specimen 


The  rivet’s  stiffness  was  calculated  using  Swift’s  displacement  compatibility  method.  In 
previous  work  [4,5]  Svrift  used  an  empirical  relationship  to  estimate  the  amount  of  load 
transfer  present  in  a  fastened  doubler.  This  method  has  also  been  shown  to  be  an  effective 
way  to  estimate  load  transfer  in  the  layered,  finite  element  code  FRANC2D/L  [6].  Even 
though  the  specimens  examined  in  this  study  are  not  of  a  “doubler-type”  configuration 
explicitly,  the  symmetry  modeling  approach  used  in  this  study  (Figure  2)  allows  the  use  of 
this  approach  as  a  reasonable  approximation. 
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The  rivet  elements  in  FRANC2D/L  are  simulated  using  an  elastic  spring  under  shear 
loading  and  each  portion  of  the  skin  and  splice  plate  is  an  idealized  bar  element.  Bar 
displacements  are  obtained  using  the  following  equation: 


PL 

AE 


where: 


P  =  bar  load 
L  =  bar  length 
A  =  bar  area 
E  =  bar  modulus 


The  rivet  stiffiiess  was  calculated  using  the  following  empirical  equation: 


f  D  D\ 

A  +  B 

— +  - 

1 _ 

where: 

F=  applied  load  (psi) 

D  =  rivet  diameter  (in) 

E  =  elastic  modulus  of  sheet  (psi) 
h  =  doubler  thickness  (in) 
ts  =  skin  thickness  (in) 

A  =  1.666  for  steel  fasteners 
B  =  0.86  for  steel  fasteners 


The  amount  of  load  transfer  present  in  a  given  specimen’s  rivet  could  then  be  calculated 
using  the  following  relationship: 


Yo  Load  Tranter  = 


Calculated  rivet  load 
Applied  remote  load 


xlOO 


This  approach  was  used  for  both  the  single  fastener  and  multiple  fastener,  in-line  specimen 
arrangements.  These  are  the  only  two  configurations  for  which  Swift’s  displacement 
compatibility  method  could  be  applied. 

For  the  multiple  fastener,  in-line  specimen,  a  “strip”  of  the  joint  was  modeled  in  the 
manner  described  above.  The  specimen  was  also  modeled  using  symmetry  conditions 
through  the  specimen’s  thickness.  The  modeled  area  of  the  specimen  is  shown  below  in 
Figure  4. 
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Area  Modeled 


Plane  of  Symmetry 


•  •  • 
•  •  • 


Plane  of  Symmetry 


Figure  4:  Area  Modeled  Using  FRANC2D/L 

A  FRANC2D/L  finite  element  analysis  was  performed  to  estimate  the  amount  of  load 
transfer  present  in  the  critical  fastener  of  the  multiple  fastener,  in-line  joint  specimen.  The 
FRANC2D/L  model  is  shown  below  in  Figure  5. 

The  results  of  the  strip  model  indicated  that  approximately  53.3%  of  the  load  was 
transferred  in  the  top  rivet  and  46.7%  of  the  load  was  transferred  from  the  bottom  rivet. 
Using  these  values  as  an  approximation,  the  top  three  fasteners  (including  the  fastener  hole 
with  the  crack)  were  estimated  to  transfer  approximately  17.7  %  of  the  load,  while  the 
bottom  fasteners  were  estimated  to  transfer  approximately  15.5%  of  the  load.  Even 
though  these  numbers  would  be  more  representative  of  an  infinitely  wide  panel  (minus 
edge  effects),  it  is  believed  that  these  percentages  would  be  sufficient  for  providing  input 
for  crack  growth  life  predictions  of  loaded  fastener  holes. 

The  multiple  fastener,  staggered  joint  specimen  arrangement  was  estimated  to  have 
approximatel\  20%  load  transfer  per  fastener  (for  life  predictions),  with  the  fasteners 
equally  sharing  the  load.  This  is  a  rather  general  approximation,  due  to  the  fact  that  there 
was  no  strict  control  on  the  finger  tight  torque  between  the  fasteners  and  variation  in  fit 
between  the  fasteners,  splice  plates  and  joint  members.  Also,  the  edge  effects  in  the 
specimens  should  have  potentially  raised  the  amount  of  load  transferred  in  the  “outer” 
fasteners  of  the  lower  row.  This  staggered  bolt  pattern  should  have  eliminated  any  load 
interaction  or  “shadowing”  between  the  two  rows  of  bolts  in  each  joined  member.  Since 
Swift’s  displacement  compatibility  method  has  only  been  shown  to  be  valid  for  in-line 
fastener  arrangements,  it  was  not  employed  here.  It  is  believed  that  this  value  (20%  load 
transfer)  is  conservative  and  should  be  sufficient  for  use  with  “loaded  hole”  crack  growth 
life  prediction  analyses  of  staggered  fastener  patterned  joint  specimens. 
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upper  Rivet  Location 


Lower  Rivet  Location 


N 


U 


--5^ 


Material:  7075-T7351  Aluminum 
E=  10.3  E6psi/in^ 
v  =  0.33 

Member  Thickness  =  0.045  in 
Splice  Plate  Thickness  =  0.09  in. 
Rivet  Stiffiiess  =  2.91  xlO^  Ib/in 


Figure  5:  Strip  Model  of  Multi  Fastener,  In-Line  Specimen 
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4. 


VERIFICATION  TESTING 


All  verification  testing  was  performed  in  the  Fatigue  and  Fracture  Test  Facility,  Bldg.  65, 
Area  B,  WPAFB,  OH.  Five  servo-hydraulic  fatigue  test  frames  were  used  to  test  titanium 
and  steel  specimens.  The  test  frames  were  operated  in  load  control  with  MTS  458  test 
controllers  at  frequencies  of  6-10  hz.  Sinusoidal  load  control  signals  were  generated  with 
MS-DOS  based  computers  running  MATE  software. 

The  test  matrix  of  40  specimens,  listed  in  Table  1,  consisted  of  center  flaw,  open  hole  and 
joined  specimens.  All  specimens  were  tested  in  lab  air.  A  more  detailed  list  including 
specimen  ID,  testing  dates  and  machine  usage  can  be  found  in  Appendix  A. 

Table  1:  Test  Matrix. 


Specimen 

Max  Load 

Stress  Ratio 

Unflawed 

Type 

(kips) 

(R) 

Specimens 

Center  Cracked 

decreasing  K 

-0.1 

2 

- 

TI-6AL-4V 

-0.5 

2 

- 

0.25  inch 

0.1 

2 

- 

thick 

0.5 

2 

- 

Center  Cracked 

decreasing  K 

-0.1 

2 

- 

4340-195ksi  Steel 

-0.5 

2 

- 

0.25  inch 

0.1 

2 

- 

thick 

0.5 

2 

- 

Comer  Cracked  Open  Hole 

16.0 

2 

- 

TI-6AL-4V 

20.0 

2 

- 

0.25  inch  thick 

24.0 

0.1 

2 

- 

Single-Fastener 

6.0 

3 

- 

Joint  7075-T73 

8.0 

■BH 

- 

3 

Multi-Fastener 

8.0 

3 

- 

In-Line  Joint  7075-T73 

11.0 

HBI 

- 

3 

Multi-Fastener 

8.0 

HBH 

3 

- 

Staggered  Joint  7075-T73 

11.0 

- 

3 

4.1.  Center  Cracked  Panel  Testing 

Titanium  (Ti-6-4  AMS  491 IG  annealed)  and  4340-195ksi  steel  panels  (3.95  x  16  x  0.25 
inch)  with  a  center  through  the  notch  (0.14  X  0.05  inch)  tested  with  a  decreasing  stress 
intensity  factor  method  in  order  to  obtain  fatigue  crack  growth  rates  of  10'*  in/cycle.  A 
saw  cut  extended  about  0.02  inch  from  the  end  of  the  notch  to  facilitate  cracking.  The 
specimens  were  fatigued  at  a  load  of  24  kips  until  cracking  began.  The  load  was 


9 


decreased  at  about  7-10%  increments  and  a  half  crack  extension  of  0.010”  until  the  crack 
arrested.  The  load  was  then  increased  at  increments  of  7-10%  until  fracture. 


4.1.1.  Ti-6-4  Center  Cracked  Panel  Testing 

Eight  titanium  specimens  (two  each  at  four  different  R-ratios)  were  tested.  Plots  of  da/dN 
versus  AK  for  R=-0.5,  -0. 1,  0. 1,  and  0.5  can  be  seen  in  Figures  6-9,  respectively.  Note; 
In  all  cases  where  R<0,  the  value  AK  is  Kn„x  since  there  is  no  physical  way  to  define  a 
negative  value  for  K. 


l.E-02 
l.E-03 
l.E-04 
l.E-05 
l.E-06 
l.E-07 
l.E-08 
l.E-09 
l.E-10 

1  10  100 


■  Ti-MT-03 
®  Ti-MT-04 


Kmax  (ksi(in)° 


Figure  6:  Crack  Growth  Rate  Plot  of  Ti-6-4  Center  Cracked  Panels  Using  the  Decreasing  K  Method 

with  R=-0.5. 
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da/dN  (in/cycle)  da/dN  (in/cycle) 


l.E-02 
l.E-03 
l.E-04 
l.E-05 
l.E-06 
l.E-07 
l.E-08 
l.E-09 
l.E-10 

1  10  100 


A 

^ - 

•An®* 

♦  ♦ 

_ 

> 

♦ 

V 

♦ 

♦  Ti-MT-02 
®  Ti-MT-10 


''igure  7: 


BCmax  (ksi(in)° 

Crack  Growth  Rate  Plot  of  Ti-6-4  Center  Cracked  Panels  Using  the  Decreasing  K 

Method  with  R=-0.1. 


l.E-02 
l.E-03 
l.E-04 
l.E-05 
l.E-06 
l.E-07 
l.E-08 
l.E-09 

1  10  100 
AK  (ksi(in)°'^) 

Figure  8:  Crack  Growth  Rate  Plot  of  Ti-6-4  Center  Cracked  Panels  Using  the  Decreasing  K 

Method  with  R=0.1. 
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da/dN  (in/cycle) 


l.E-02 
l.E-03 
l.E-04 
l.E-05 
l.E-06 
l.E-07 
l.E-08 
l,E-09 

1  10  100 


e 

o 

o® 

m 

o 

■  Ti-MT-08 
®  Ti-MT-09 


AK  (ksi(in)°'^) 

Figure  9:  Crack  Growth  Rate  Plot  of  Ti-6-4  Center  Cracked  Panels  Using  the  Decreasing  K 

Method  with  R=0.5. 


A  composite  of  all  the  Ti-6-4  data  is  shown  in  Figure  10,  Each  series  of  data  includes  the 
two  specimens  represented  in  Figures  6-9.  These  data  were  used  to  modify  AFGROW’s 
material  data  file. 
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l.E-02 
l.E-03 
l.E-04 
■R  l.E-05 

o 

i  l.E-06 
1  l.E-07 

T3 

l.E-08 
l.E-09 
l.E-10 

1  10  100 

AK  (ksi(in)°'^) 

Figure  10:  Composite  of  All  Ti-6-4  Center  Cracked  Panel  Tests.  Each  Series  of  R-ratio  Data 

Represent  Two  Specimens. 

4.1.2.  4340-195ksi  Steel  Center  Cracked  Panel  Testing 

Eight  steel  specimens  (two  each  at  four  different  R-ratios)  were  tested.  Plots  of  da/dN 
versus  AK  for  R=-0.5,  -0.1,  0.1,  and  0.5  can  be  seen  in  Figures  11-14,  respectively. 
Note:  In  all  cases  where  R<0,  the  value  AK  is  Kmax  since  there  is  no  physical  way  to 
define  a  negative  value  for  K. 
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da/dN  (in/cycle)  da/dN  (in/cycle) 


l.E-02 
l.E-03 
l.E-04 
l.E-05 
l.E-06 
l.E-07 
l.E-08 
l.E-09 

1  10  100  1000 
Kmax  (ksiCin)”'^) 

Figure  11:  Crack  Growth  Rate  Plot  of  4340-195ksi  Steel  Center  Cracked  Panels  Using  the 

Decreasing  K  Method  with  R— 0.5. 

l.E-02 
l.E-03 
l.E-04 
l.E-05 
l.E-06 
l.E-07 
l.E-08 
l.E-09 

1  10  100  1000 
K„,ax  (ksi(in)°'’) 

Figure  12:  Crack  Growth  Rate  Plot  of  4340-195ksi  Steel  Center  Cracked  Panels  Using  the 

Decreasing  K  Method  with  R=-0.1. 
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da/dN  (in/cycle)  da/dN  (in/cycle) 


l.E-02 
l.E-03 
l.E-04 
l.E-05 
l.E-06 
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Figure  13:  Crack  Growth  Rate  Plot  of  4340-195ksi  Steel  Center  Cracked  Panels  Using  the 

Decreasing  K  Method  with  R=0.1. 
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Figure  14:  Crack  Growth  Rate  Plot  of  4340-195ksi  steel  Center  Cracked  Panels  Using  the 

Decreasing  K  Method  with  R=0,5, 
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The  data  for  St-MT-07  in  Figure  14  are  nonlinear  because  the  load  was  incrementally 
increased  greater  than  10%  after  crack  arrest.  A  composite  of  all  the  steel  data  is  shown  in 
Figure  15.  Each  series  of  data  includes  the  two  specimens  represented  in  Figures  11  -  14. 
There  appears  to  be  two  rather  distinct  curves  in  Figure  15.  The  curve  which  lies  in  the 
intersection  of  AK=10ksiVin  and  da/dN=lE-5  is  very  atypical  of  steel  data  [7].  The  data 
were  also  compared  to  test  data  for  4340-180ksi  steel  to  verify  that  the  test  points  taken 
while  the  load  increase  was  high  were  unreasonable.  These  data  were  removed  from  the 
plots  and  the  result  is  given  in  Figure  16.  These  data  were  used  to  modify  AFGROW’s 
material  data  file. 

l.E-02 

l.E-03 

^  l.E-04 

^  l.E-05 

l.E-06 

1  l.E-07 
•a 

l.E-08 
l.E-09 
l.E-10 

1  10  100  1000 

AK  (ksi(in)° 

Figure  15:  Composite  of  All  4340-195ksi  Steel  Center  Cracked  Panel  Tests.  Each  Series  of  R-ratio 

Data  Represents  Two  Specimens. 
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z  l.E-06 
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l.E-08 


l.E-09 

1  10  100  1000 


♦R=-0.5 

®R=-0.1 

^R=0.1 

XR=0.5 


AK  (ksi(in)°'^) 

Figure  16:  Same  Graph  as  Figure  15  with  Unreasonable  Data  Removed. 


4.2.  Open  Hole  Ti-6-4  Panel  Testing 

Six  open  hole  titanium  specimens  were  tested.  A  0.02  to  0.03  inch  (length  and  depth) 
comer  notch  was  cut  on  one  side  of  the  0.1875  inch  diameter  hole.  The  specimen  was 
precracked  at  a  load  of  33  kips  until  the  crack  grew  0.06  inch  beyond  the  hole.  The  hole 
was  then  redrilled  to  0.25  inch  diameter  and  the  specimen  was  cyclically  loaded  at  the 
same  stress  level  until  the  crack  length  was  0.03  to  0.05  inch  beyond  the  new  hole 
diameter.  The  specimens  were  cyclically  load  shed  at  7-10%  load  increments  and  0.1  inch 
crack  increments  until  the  testing  load  was  reached.  Two  specimens  each  at  loads  of  16, 
18,  and  20  kips,  R=0.1  were  tested.  The  crack  length  was  recorded  at  0.010  inch 
increments  to  determine  when  the  comer  crack  transitioned  to  a  through  crack.  Once  the 
crack  became  a  through  crack  the  crack  length  was  recorded  at  0. 1  inch  increments  until 
failure. 


The  crack  growth  life  was  compared  to  AFGROW  predictions.  AFGROW  v3.73  was 
used  with  the  updated  material  data  file  including  data  generated  for  this  report.  The 
updated  material  data  file  was  used  to  make  these  predictions.  A  comparison  between  the 
test  data  (from  Figure  8)  and  the  fit  used  in  the  AFGROW  database  for  Ti-6-4  AMS 
491 IG  annealed  is  given  in  Figure  16.  Figures  17-22  show  the  crack  growth  data  for  each 
specimen  with  the  AFGROW  prediction.  The  transition  from  a  comer  crack  to  a  through 
crack  is  also  noted  on  each  graph.  The  initial  crack  length  was  known  at  the  beginmng  of 
each  test,  but  the  crack  depth  was  not  known.  The  crack  depth  was  estimated  using 
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AFGROW.  A  0.025  inch  quarter-circle  corner  crack  was  grown  from  a  0.1875  inch  hole 
until  crack  equaled  the  original  crack  length  after  the  hole  was  redrilled.  Due  to  re-drilling, 
the  aspect  ratio  of  the  crack  changed.  The  crack  depth  was  estimated  by  multiplying  the 
depth  of  the  crack  on  the  0.1875  inch  hole  by  0.886  and  inputting  the  result  as  the  crack 
depth  on  the  0.25  inch  hole.  The  crack  was  then  grown  to  emulate  the  load  shedding  that 
occurred  before  actual  test  data  were  taken.  The  initial  crack  sizes  input  to  AFGROW  are 
located  in  Table  2.  AFGROW  was  then  restarted  and  the  results  are  plotted  in  Figures  17 
-  22  along  with  the  data  for  each  specimen.  In  some  cases,  a  second  crack  formed  on  the 
opposite  side  of  the  hole  during  a  test.  The  effect  of  this  cannot  currently  be  predicted 
using  AFGROW.  However,  it  will  cause  acceleration  in  growth  of  the  primary  crack. 


Crack  growth-rate  data 


Figure  17:  Test  Data  Shown  in  Figure  8  Plotted  with  AFGROW  Database  Plot  for  Various  R-ratios. 


Table  2:  Crack  Sizes  Input  to  AFGROW  for  Life  Prediction  Shown  in  Figures  17  -  22. 


Specimen 

Crack  length,  c 

Crack  depth,  a 

Ti-OH-01 

0.087 

0.128 

Ti-OH-02 

0.103 

0.146 

Ti-OH-03 

0.047 

0.080 

Ti-OH-04 

0.049 

0.082 

Ti-OH-05 

0.055 

0.094 

Ti-OH-06 

0.081 

0.119 
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♦  Ti-OH-1  Data 
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W  Breakthrough _ 


Figure  18:  Constant  Amplitude  (16  kips)  Testing  Data  of  Ti-OH-1  and  AFGROW  Prediction. 


0  25,000  50,000  75,000  100,000  125,000  150,000  175,000 

Cycles 


Figure  19:  Constant  Amplitude  (16  kips)  Testing  Data  of  Ti-OH-2  and  AFGROW  Prediction. 
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♦  Ti.0H.3Data 
AFGROW  Prediction 

*  Breakthrough _ 


igure  20: 


Constant  Amplitude  (18  kips)  Testing  Data  of  Ti-OH-3  and  AFGROW  Prediction. 


Cycles 


Figure  21:  Constant  Amplitude  (ISkips)  Testing  Data  of  Ti-OH-4  and  AFGROW  Prediction. 
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♦  Ti-OH-5  Data 
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*  Breakthrough _ 


i^'igure  22:  Constant  Amplitude  (20  kips)  Testing  Data  of  Ti-OH-5  and  AFGROW  Prediction. 


Cycles 


Figure  23:  Constant  Amplitude  (20  kips)  Testing  Data  of  Ti-OH-6  and  AFGROW  Prediction. 
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4.3.  Fastener  Joint  Specimens 


Flawed  and  unflawed  fastener  joint  specimens,  seen  in  Figure  1,  were  fatigue  tested  under 
constant  amplitude  conditions  until  failure.  The  flawed  specimens  were  precracked  in  the 
same  manner  as  the  open-hole  titanium  specimens.  The  panels  with  a  single  0.1875  inch 
open  hole  and  an  initial  comer  saw  cut  of  0.02-0.03  inch  length  and  depth  were 
precracked  to  a  length  of  0.07  inch.  The  panels  were  then  cut  in  half  to  a  length  of  8 
inches  and  the  precracked  hole  was  then  redrilled  to  0.25  inch  diameter  (leaving  a  0.04 
inch  precrack).  A  schematic  of  a  specimen  is  shown  in  Figure  23.  The  panels  were  then 
riveted  together  as  shown  in  Figure  1 .  The  cracked  hole  was  in  the  top  middle  hole  of  the 
in-line  joint  specimen  and  in  one  of  the  top  holes  with  the  crack  growing  towards  the  other 
rivet  in  the  staggered  joint  specimen.  The  test  matrix  and  lifetimes  are  shown  in 
Appendix. 
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0.1875”  hole 
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0.25”  holes 
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o  o  o 

5” 
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Figure  24:  Schematic  of  the  Joint  Specimens  After  Precracking  but  Before  Riveting. 


4.4.  AFGROW  Predictions  of  Joint  Specimens 

Using  initial  conditions,  AFGROW  was  used  to  predict  the  lifetimes  of  the  specimens.  The 
test  results  were  compared  to  AFGROW  predictions.  The  AFGROW  predictions  were 
handled  in  the  same  manner  as  the  open  hole  specimens.  An  initial  saw  cut  was  modeled 
as  a  quarter  circle  comer  crack  in  a  0.1875  inch  diameter  hole.  The  crack  was  then  grown 
to  the  known  length  before  the  hole  was  drilled  to  0.25  inch.  Then  the  test  prediction 
began.  In  order  to  predict  the  life  of  a  joint  specimen  with  AFGROW,  the  specimen 
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configuration  must  first  be  modeled.  The  single  joint  configuration  was  easy  to  model 
with  AFGROW.  The  comer  crack  at  open  hole  configuration  was  chosen  with  a  load 
transfer  of  100%.  (100%  of  the  load  was  introduced  through  the  fastener.) 


Modeling  the  load  transfer  coefficient  for  the  staggered  and  in-line  joint  specimens  is  more 
complicated  however.  To  further  illustrate  the  load  transfer,  consider  Figure  24.  The 
stresses  at  the  edge  of  the  holes  at  cut  A  result  from  bearing  stresses  from  the  loaded 
fasteners  in  cut  A  as  well  as  bypass  stresses  from  the  loaded  fasteners  in  cut  B.  The  stress 
at  the  edge  of  the  holes  at  cut  B,  however,  is  a  result  of  bearing  stresses  only  from  the 
loaded  fasteners  in  cut  B.  The  equation  for  determining  the  load  transfer  coefficient  is 


given  by. 


3/> +2P, 


or 


i.r.B  = 


3P. 


t  t  t  t  t  , 


A 

_ _ jr\ _ px _ 

B 

W  Vy'  Fa 

Pb  Pb  Pb 

Figure  25:  Schematic  of  Staggered  Joint  Specimen  and  Loads  Used  to  Calculate  Load  Transfer 

Coefficient. 

The  load  transfer  per  fastener  was  discussed  in  Section  1,  and  for  the  staggered  joint 
specimens,  the  load  transfer  per  fastener  was  determined  to  be  20%.  Since  only  two  of 
the  five  rivets  of  the  staggered  joint  specimens  are  in  the  top,  the  load  transfer  coefficient 
is  40%.  The  staggered  joint  specimens  were  modeled  as  an  open  hole  with  a  double 
comer  crack.  Although  only  one  of  the  holes  was  precracked,  examination  of  the  fracture 
surface  showed  that  a  crack  appeared  to  grow  from  the  other  hole  at  an  equal  rate  and 
coalesced,  and  the  two  cracks  coalesced  midway  between  the  holes.  For  this  reason  the 
specimen  was  assumed  to  be  fractured  when  the  crack  length  was  0.375inch.  The  width  of 
the  specimens  was  assumed  to  be  1.75  inch  since  that  is  half  the  distance  to  the  nearest 
hole.  Although  the  width  of  the  specimen  was  3.95  inches,  a  width  of  1.75  inches  was 
used  because  the  nearest  hole  was  assumed  to  be  an  edge.  This  engmeenng 
approximation  is  reasonable  because  as  the  crack  approaches  the  hole,  it  will  speed  up  as 
if  it  is  approaching  an  edge,  and  once  the  crack  reaches  the  adjacent  hole,  the  rest  of  the 
life  is  short  compared  to  the  time  it  took  to  reach  the  hole. 


23 


For  the  in-line  joint  specimens,  the  specimens  were  modeled  as  a  single  hole  with  a  comer 
crack.  The  top  row  of  in-line  rivets  shadow  the  bottom  row,  which  results  in  the  top  row 
transferring  53.3%  of  the  load,  as  discussed  in  the  finite  element  section.  The  adjacent 
hole  was  0.875  inch  fi'om  the  cracked  hole,  and  the  width  of  the  modeled  specimens  was 
twice  that  distance,  1.75  inches. 

The  predicted  fatigue  lives  are  compared  to  the  actual  lifetimes  and  the  results  are  shown 
in  Table  3.  The  single,  in-line,  and  staggered  joint  specimens  had  an  average  error  of 
18%,  4%,  and  19%.  These  results  are  promising  considering  that  AFGROW  currently 
does  not  have  the  capability  to  model  multiple  hole  specimens. 


Table  3:  Flawed  Riveted  Specimen  Lifetimes  and  Information  Used  for  AFGROW  Predictions 

Where  c  is  Measured  and  a  is  Calculated. 


Specimen 

■Esm 

c  (in.) 

a  (in.) 

Lifetime 

Prediction 

4-04 

6 

0.036 

0.0779 

416,632 

274,073 

4-05 

6 

0.039 

0.0814 

302,224 

268,155 

4-06 

■B9 

6 

236,673 

256,000 

■EBEEH 

in-line 

8 

0.040 

■Mils™ 

151,008 

146,461 

in-line 

8 

0.039 

WbtiIMM 

139,660 

146,725 

8 

0.039 

0.0684 

54,945 

75,000 

5B-08 

staggered 

8 

0.042 

0.0710 

73,540 

74,000 
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5. 


AFGROW  ENHANCEMENTS 


Various  features,  including  residual  stress  effects,  multi-dimensional  crack  growth  rate 
data,  user  input  beta  factor  table,  and  real-time  plotting  have  been  added  to  AFGROW. 

The  multi-dimensional  crack  growth  rate  data  have  been  added  to  predict  the  crack 
growth  life  of  a  structure  which  may  experience  a  variety  of  environments  during  its 
lifetime.  Crack  growth  rate  data  can  be  stored  and  retrieved  for  a  given  material  imder 
different  enviromnental  conditions.  The  data  can  be  accessed  according  to  their 
environmental  condition  as  a  function  of  crack  position.  The  enhancement  has  the 
following  specifics:  1)  user  can  input  modified  crack  growth  rated  data  for  a  given 
material  under  different  environmental  conditions,  2)  user  can  save  crack  growth  data  to  a 
file  or  read  crack  growth  data  from  a  file,  3)  user-input  crack  growth  rate  data  can  be  used 
with  both  tabular  material  data  and  Walker  equation  options,  4)  user  can  specify 
environmental  condition  and  region  around  a  crack  as  a  function  of  distance  “c”  (surface 
direction)  and  “a”  (thickness  direction,  where  applicable)  ahead  of  the  crack  front(s),  5) 
AFGROW  will  transition  from  one  environmental  condition  to  another  by  a  polynomial 
function  (to  third  order),  as  the  crack  grows.  The  environment  dialog  box  and  specimen 
model  with  tw^o  environments  can  be  seen  in  Figures  25  and  26,  respectively. 


iocafioti  of  Environment 
Locations  are  measured  from  center  of  crack 

Oir 


Figure  26:  Environment  Dialog  Box. 
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Double  corner  crack  at  hole 


d  =  0.25 


w  =  4.00 


Environments: 

B  z5 
B  pure 


Figure  27:  Model  With  Two  Environments  Applied. 

The  Beta  correction  table  allows  the  user  to  enhance  the  standard  solutions  currently  used 
to  calculate  the  crack  driving  parameter,  K.  For  example,  consider  a  crack  emanating 
from  a  notch  with  a  large  stress  concentration.  First  to  model  this  situation,  a  standard 
crack  configuration  is  chosen,  then  the  Beta  correction  table  is  used  to  modify  the  crack 
driving  force  as  a  function  of  distance  from  the  notch.  More  specifically,  the  input  Beta 
Factor  table  will  operate  as  such:  1)  AFGROW  default  will  have  the  Beta  Factor 
modification  turned  off,  2)  the  user  input  Beta  Factor  table  could  be  input  manually,  or 
read  from  a  file,  3)  the  user-input  Beta  Factor  table  function  can  be  turned  off/on  at  any 
time.  In  addition,  the  graphical  user  interface  will  be  written  so  that  the  appropriate 
portions  of  the  code  will  be  disabled  when  the  capability  is  used,  and  that  existing 
methods  remain  in  the  code  after  this  modification  is  made.  The  Beta  Factors  dialog  box 
can  be  seen  in  Figure  27. 
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BHA  fACTOfiS 

Nonmtltjod  Stress  thstnlMJtftxis 
///  ttie  fvUuwifnj  furnjdt 

Use  itff  to  ?‘i  sets  of  for  ffXTn^sfog  'p '  fits tof ices 

At/  a  Stress(x)  iifid  Strvss{y)  trivc  Been  set  to  TO 
Sabseqoont  vetoes  ore  norms^/fxi  yert  Stress  otz-Q 


Natobcr  of  Sets:  t 

Set  /Coord  Stress(X)  StressfY) 

0  (10  r.Q  to 


Figure  28:  Beta  Factor  Table. 


Real  time  plotting  and  output  file  creation  have  been  added  to  AFGROW.  The  real  time 
plotting  capability  plots  the  crack  length  versus  cycles  vyhile  AFGROW  is  running.  Up  to 
eight  different  conditions  can  be  overlaid  on  the  same  plot.  The  output  file  contains  data 
columns  of  crack  length,  cycles,  stress  intensity,  beta,  stress  ratio,  stress,  and  crack 
growth  rate.  The  plot  window  can  be  seen  in  Figure  28. 
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Figure  29:  Plot  Window  Showing  Crack  Length  Versus  Cycles. 
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7.  APPENDIX 


7.1  Detailed  Specimen  Information 


^-Pecrtaaing'K Testae  _ 


Spec  ID  »\ 

Mach  W  f 

t  1 

oad  (kips  ’ 

R 

Begin 

End  ! 

TI-MT-02 ! 

13  ’  3.9505  1 

0.26040 

24.9  i 

-0.1 

07-19-95 

Overload 

Ti-MT-03 

12  :  3.9505  : 

0.25820 

24  i 

-0.5 

11/3/95 

4/9/96 

Ti-MT-04 

4  i  3.9500  i 

0.25240  ; 

24 

-0.5  ; 

11-15-95; 

3/27/96 

Ti-MT-06  i 

12  3.9500  I 

0.26480 

26 

0.1 

07-24-95  !  11-03-95  i 

Ti-MT-07  ' 

13  3.9500  j 

0.25910 

24 

0,1 

10-25-95; 

3/26/96  : 

Ti-MT-Oe : 

10  ,  3.9500 

0.25680  ; 

24 

0.5 

11-20-95; 

4/12/^ 

Accidental  Crush 

Tt-MT-09 

1  '  3.9500  1 

0.25710 

24 

0.5 

11-15-95  i 

6/13/96 

Ti-MT-10 

13  3.9500  , 

0.25900 

24  ! 

^.1  i 

3/27/96'! 

5/17/96  : 

I 

1 

1 

|i 

1 

1. 

Decreasing  K  testing; 

: 

Spec  ID  0 

Mach  W 

t 

oad  (kips 

R 

Begin 

End 

St-MT-01 

iot  3.945! 

0.213 

30 

-0.5 

4/24/96 

5/24/96 

cf  not  straight 

St-MT-02 

ioi  3.95 

0.212 

30 

4D.5 

5/30/96 

7/16/96 

St-MT-03 

12:  3.9485 

0.211 

30  1 

-0.1 

4/25/96 

6/11/96 

St-MT-04 

12  3.94 

0.225 

30 

-0.1 

6/12/96 

7/2/96 

St-MT-05 

4i  3.9465 

0.214 

30 

0.1 

5/2/96 

5/30/96 

not  grow  -grow  ! 

St-MT4D6 

4!  3.95^ 

0.226' 

30 

0.1 

5/30/96 

8/5/96 

! 

St-MT-07 

13  3.95 

0.208 

30 

0.5 

6/7/96 

7/25/96  i 

St-MT-08 

12'  3.952 

0.235 

30 

0.5 

7/2/96 

8/20/96  ' 

LCobitalnt  L^Testing 

Spec  ID  0 

Mach  I  W 

t 

oad  (kips 

R 

Cycles 

Begin 

End 

Ti-OH-01 

4  !  3.956 

0.251 

16 

0.1 

233,650 

4/17/96 

4/19/96 

Ti-OH-02 

4  '  3.946  “ 

0.257 

16 

0.1 

141,226 

4/22/96 

4/23/96 

Ti-OH-03 

10  3.9615 

0.24 

20 

0.1 

67,450 

4/18/96 

4/18/96 

:  Ti-OH-04 

10  i  3.9465 

0.24 

20 

0.1 

71,750 

4/19/96 

4/19/96 

I  Ti-OH-05 

10  3.947 

0.249 

24 

0.1 

38.577 

4/22/96 

i  4/22/96 

I  Ti-OH-06 

10  ;  3.949 

0.245 

24 

0.1 

33,320 

4/23/96 

1  4/23/96 

1 

A  . 

[Spec  IDff 

Mach  W 

t 

PreCrack 

Stress  (ksi) 

Cydes 

1  4-04 

13  1.0055 

0.250 

0.036 

6.0 

[  416,632 

1  4-05 

13  i  1.0050 

0.250 

!  0.039 

6.0 

302,224 

approx  (Power  Out) 

1  4-06 

13  1.0050 

0.250 

i  0.048 

6.0 

;  236,673 

. . .  --S-iir: 

[Spec  ID  #1  Mach  Wit 

PreCrack 

Stress  (ksi) 

1  Cycles 

!  4-01 

12  '  1.0049 

;  0.250 

** 

8.0 

!  400,211 

'  4-02 

1  12  1.0050 

0.250 

** 

:  8.0 

^  244,675 

,  4-03 

i  12  :  1.0050 

:  0.250 

** 

8.0 

235,616 

iapprox  (Power  Out) 

Tbwed  Mudi  Fastener  In  Line  J<^t 

. 

......... 

. 

Spec  ID#!  Mach  W 

i  t 

■  PreCrack  Stress  (ksi) 

Cycles 

' 

.  5A-07 

;  12  3.9380 

1  0.090 

0.040 

8.0 

151,008 

After  Bolts  Loosened  (52,253) 

;  5A-08 

12  ;  3.9380 

0.090 

0.039 

8.0 

139,660  i  Bolls  Loose  1 

. . .  ..........,.,,5.,, 

Uhflawed  Mutli  Fi^ner  In  line  Joint 

■  '  ■ ■.'  =  ■  ■' . :.,"V 

[Spec  ID#  Mach  '  W 

t 

PreCrack  |  Stress  (ksi) 

1  Cydes 

i 

'  SA-Ol 

i  12  ;  3.9500 

1  0.090 

** 

11.0 

,  524,611 

1  Bolts  Tight 

,  5A-02 

1  12  i  3.9485 

’  0.090 

*• 

11.0 

!  98,793 

Bolts  Loose  i 

1  5A-03' 

12  :  3.9480 

,  0.090 

** 

[  11.0 

,  77,676 

[After  Bolts  Loosened  (223,627) 

'  MiM  FasttMo*  Stagi^d  Jpim _ _ ;  _  _ _ _ 

[  Spec  t  !  PreCrack  I  Stress  (ksi)  i  Cycles 


1  spec  11/ #  ;via(;n 

5B-07  ;  13 

3.9390  ! 

0.090 

0.039 

8.0 

57,945  1  After  Bolts  Loosened  (52,060) 

;  5B-08  13 

3.9350 

0.090  i 

0.042 

8.0 

73,540  Bolts  Loose 

iJbfbwt^  MvtU  Fasten^  Steggei^^^ 

[Spec  ID  #i  Mach 

W 

. 

PreCrack 

I  Stress  (icsij  i 

Cycles  : 

I  5B-01  1 

:  3.9490 

0.090 

** 

11.0 

45,922  'Bolts  Loose 

;  5B-02  i 

3.9490 

0.090  ; 

** 

11.0 

54,327  Bolts  Loose 

'  5B-03  : 

3.9490 

0.090  ' 

** 

1  11.0 

;  58,115  1  Bolls  Loose 

30 


